beh(1).M(1).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9719.lvd');
beh(1).M(2).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9720.lvd');
beh(1).M(3).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9721.lvd');
beh(2).M(1).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9725.lvd');
beh(2).M(2).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9726.lvd');
beh(2).M(3).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9727.lvd');
beh(3).M(1).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9728.lvd');
beh(3).M(2).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9728.lvd');
beh(3).M(3).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9730.lvd');
beh(4).M(1).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9749.lvd');
beh(4).M(2).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9753.lvd');
beh(4).M(3).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9754.lvd');
beh(5).M(1).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9755.lvd');
beh(5).M(2).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9757.lvd');
beh(5).M(3).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9756.lvd');
beh(6).M(1).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9782.lvd');
beh(6).M(2).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9783.lvd');
beh(6).M(3).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9793.lvd');
beh(7).M(1).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9802.lvd');
beh(7).M(2).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9803.lvd');
beh(7).M(3).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9804.lvd');
beh(8).M(1).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9805.lvd');
beh(8).M(2).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9806.lvd');
beh(8).M(3).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9807.lvd');
beh(9).M(1).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9808.lvd');
beh(9).M(2).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9809.lvd');
beh(9).M(3).data=load_lvd('D:\tmpTest\Behavior M086 - 089\S1-T9810.lvd');

for lnd = 1:length(beh)
    for knd = 1:length(beh(lnd).M)
        for ind = 1:length(beh(lnd).M(knd).data(10,:))
            if beh(lnd).M(knd).data(10,ind)>2.5
               beh(lnd).M(knd).angle(ind)=(beh(lnd).M(knd).data(10,ind)-2.5)/2.5*180;
            end
            if beh(lnd).M(knd).data(10,ind)<2.5
               beh(lnd).M(knd).angle(ind)=(abs(beh(lnd).M(knd).data(10,ind)-2.5))/2.5*(-180);
            end
        end
    end
end

figure,
id = 0;
for ld = 1:length(beh)
    for kd = 1:length(beh(ld).M)
        id = id+1;
        subplot(length(beh),length(beh(ld).M),id), hist((beh(ld).M(kd).angle),100)
        xlim([-185 185])
        ylim([0 400000])
    end
end

for lnd = 1:length(beh)
    for knd = 1:length(beh(lnd).M)
        beh(lnd).M(knd).polar = beh(lnd).M(knd).angle * pi/180;
    end
end

for lnd = 1:length(beh)
    figure, rose(beh(lnd).M(3).polar,1000)
end


h = [subplot(231) subplot(232) subplot(233) subplot(234) subplot(235) subplot(236)] 
axis([h], xlim)
